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We consider electronic transport through break-junctions 
bridged by a single molecule in the Kondo regime. We 
describe the system by a two-channel Anderson model. 
We take the tunneling matrix elements to depend on 
the position of the molecule. It is shown, that if the 
modulation of the tunneling by displacement is large, 
the potential confining the molecule to the central po- 
sition between the leads is softened and the position 
of the molecule is increasingly susceptible to external 
perturbations that break the inversion symmetry. In this 
regime, the molecule is attracted to one of the leads 
and as a consequence the conductance is small. We ar- 
gue on semi-classical grounds why the softening occurs 
and corroborate our findings by numerical examples ob- 
tained by Wilson's numerical renormalization group and 
Schonhammer-Gunnarsson's variational method. 
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Sketch of the model device. The overlap integrals between the 
molecular orbital and the leads (Vl,Vr for left, right lead, re- 
spectively) are modulated by the position of the molecule. The 
energy of the molecular level e is determined by the gate volt- 
age. 
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1 Introduction and model In recent years consider- 
able advance has been achieved in manufacturing meso- 
scopic systems, which due to their tunability with external 
electrodes provide a playground for research in the corre- 
lated electron systems. For example, the Kondo effect - 
generic name for phenomena related to increased scatter- 
ing off impurities with internal degrees of freedom - was 
observed in measurements of electron transport through 
quantum dots [1|, atoms, and molecules |2 3 4 5 6][7][8]. 



frequencies of the molecular internal modes it was shown 
E|, that in some cases also the oscillations of the molecule 
with respect to the leads have to be taken into account. 

In this work we concentrate on such a case. The tun- 
neling is generally dependent on the overlap between the 
wave-functions. This motivates us to investigate the effects 
of the modulation of the tunneling matrix elements by molec- 
ular oscillations. 



Unlike in experiments with quantum dots, the transport 
through break junctions is strongly affected by the molecu- 
lar vibrational modes, because the frequency of the oscilla- 
tions is of comparable magnitude than other energy scales, 
such as Coulomb repulsion. For example, the side-peaks in 
the non-linear conductance [4 5 6 1 were observed indicat- 
ing the transfer of energy from the oscillations to the elec- 
tron current. By comparing the observed frequencies to the 



More specifically, we describe the break junction by a 
model consisting of two metallic leads (half-filled bands 
of non-interacting electrons). The leads are bridged by a 
single molecule which we assume is confined harmoni- 
cally to the center between the leads. The position of the 
molecule x determines the tunneling matrix elements. As- 
suming the leads are identical and the displacement of the 
molecule from the center of inversion x is small [9|, the 
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tunneling matrix elements towards left and right leads read 
V L , R (x) = (lT9x)V. 

The Hamiltonian consists of several parts 



H = Hr 



H mo \ + H, 



vib 



H' . 



The first two terms describe the isolated leads 



E 



tkC \aa Ck ' 



(1) 



(2) 



a = L,R is the lead index, a =T, j is the spin index 
and k the wave-vector index; c* (as well as S to be intro- 
duced below) denote the fermion creation operators. The 
precise form of the band dispersion is not important for 
the results presented here, provided the density of states 
is smooth and symmetric around the Fermi surface; never- 
theless for definiteness we note that in our calculations we 
used flat (NRG simulations) and tight-binding bands (sim- 
ulations with Schonhammer-Gunnarsson wave functions). 

The second term describes the molecular orbital (we 
assume a single molecular level participates actively in the 
transport), 



(3) 



where n a = d\d a counts the number of electrons occu- 
pying the orbital, e the energy of the orbital relative to the 
chemical potential of the leads, U is the Coulomb repul- 
sion. We concentrate on the particle-hole symmetric case 



-U /2. i?vib = da) a describes the phonon mode ( 



is the phonon creation operator). 

The coupling between the leads, the molecular orbital 
and also phonons is described by 



h' = vJ2 [(i-C-sz)cLA- 



(4) 



+ (1 + C + gzWkRcrda + h ' C 



Here we introduced a constant £, which we will use to test 
for the influence of the breaking of the inversion symmetry; 
for C = the Hamiltonian is symmetric with respect to 
operation x — > —x,L<-> R; finite £ breaks this symmetry. 

In the model with no coupling to phonons (g = 0) 
for £ — only the even combination of the operators, 
i.e. c e = (cl + cr)/V2 (other indeces are suppressed) 
in the leads is coupled to the molecular orbital. The life 
time of electrons on the orbital is finite due to the tun- 
neling to the leads; the hybridization r (inverse life-time) 
for the flat band reads r = 2npV 2 , where p is the den- 
sity of states in each of the two leads (for flat band of 
half-width D, p = 1/(2D). The odd linear combinations 
c o = ( c l — c Ft,)/V2 are decoupled. For g = and £ > 
still a particular linear combination of the states in the leads 
is decoupled and the system can be described by the single- 
channel model. 



For finite g ^ this no longer holds and the molecule 
is coupled to both conduction channels. Rewriting the cou- 
pling term in the even-odd basis, 



H' = V2VJ2 [4 e< A + (C + 9x)c\ oa d 

ka 



h. 



(5) 



makes it manifest that we are dealing with a two-channel 
Anderson model. No linear combination of the conduction 
electrons can be integrated out because the coupling to the 
odd channel is mediated by phonons. This occurs because 
the modulation of tunneling is antisymmetric with respect 
to inversion; if the modulation of tunneling is symmetric 
I10 H] the orbital is still coupled only to the even channel. 

It is difficult to access the low temperature (Kondo) 
regime because of the presence of exponentially small en- 
ergy scale Tr oc cxp [— l/{pj)] ( J ~ V 2 /U is the mag- 
nitude of the exchange coupling). Only few methods [12] 
reproduce the increase of conductance towards the unitary 
limit Go = 2e 2 /h for temperatures smaller than Tk accu- 
rately. On the other hand, in the high-temperature regime 
(relevant for nanoelectromechanical systems IU3lll4|[T5l ) it 
is adequate to ignore the Kondo correlations and take only 
lowest orders in tunneling into account. 

For the work on this model in the Kondo regime, which 
was originally stimulated by Ref. lfT6l . we refer the reader 
to Refs. ||17||18|[T9l l9l. Another very recent paper discusses 
a similar model as an example of two-level system j20| . 

In the following section we first demonstrate that the 
confining potential is weakened by the electron-phonon cou- 
pling and can even be driven to a form of the double-well. 
Then we give numerical examples on how the emergence 
of the double-well potential affects static and dynamic prop- 
erties of the molecule. Finally, we show that when the in- 
version symmetry is not perfect (£ ^ 0), the electron- 
phonon coupling will drive the molecule away from the 
central position. As a consequence, in this regime the con- 
ductance can be significantly suppressed. 

2 The emergence of a double well potential It is 

easy to demonstrate that a double well effective oscilla- 
tor potential may form under the influence of a sufficient 
electron-phonon coupling g > g<j, where we define the de- 
limiting value of the coupling constant. For U = we can 
make a simple estimate on the form of the effective oscil- 
lator potential by a substitution a, ~ cv — > x/2, where x is 
a real-valued constant. In the wide-band limit (r/ D small) 
the energy gain due to the hybridization is 11211 AEh y \, = 
—2/irr log D /P, where we use the effective displacement- 
dependent hybridization r(x) = r(l + g 2 x 2 ). The elastic 
energy cost is AE C \ = fix 2 /4, hence in this semi-classical 
approximation the dependence of energy AE§c — AE e \ + 
AEhy\, on x can be written in a closed form 

AE sc (x) = f2x 2 /A - (2/7r)r(a;) Iog{D/[f(i)]}. (6) 

The prefactor of the a; 2 -term in the small-x expansion is 
equal to (2/1 - {(2/ir)g 2 r [log(D/T) - 1]}, hence for 
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Figure 1 Numerical results for effective potential: ( = 
(full lines), £ = 0.01 (dashed lines). Semi-classical esti- 
mate is also shown (dotted). Here and in subsequent fig- 
ures we use U = 0.3, r = 0.02, Q = 0.2. We use the band 
half-width D as the energy unit. The curves are shifted ver- 
tically so that the values at x = match. 



increasing g the elastic potential is softened and a double 
well effective potential emerges for 



.9 > 9d 



sr [io g (£>/r) - i] ' 



(7) 



We plot AE sc (x) for g = 0.8 < g d and g = 0.9 > g d in 
Fig- 03 (thin, dotted lines). 

We estimate the effective potential also numerically us- 
ing a variational method based on the Schonhammer-Gunnar- 
sson I22j [23l wave fuction (the details of our implementa- 
tion are given in our previous work 124II25|[TTI ). Briefly, 
the idea is to find an auxiliary non-interacting Hamiltonian 
H [of the same form as H in Eq. (Q]), but for g = 0, U = 
and renormalized parameters Vl, Vr, e], which minimizes 
the variational ground state energy E = (<I r \H\ [ I r ). The 
variational function \P is expressed in the basis of projec- 
tion operators Pi acting on the Hartree-Fock ground state 
|<?b) (which includes the phonon vacuum) of the auxiliary 
Hamiltonian H, 



W) = Y,i ) ™((J) n p l \v ). 



(8) 



To obtain the ground state, we minimize energy with re- 
spect to all the parameters of H. On the other hand, by re- 
stricting the minimization to a particular subspace (for ex- 
ample, by fixing the ratio Vl /Vr = r) we obtain the vari- 
ational wave-function $f r for which the expectation value 
(x) r is a function of r. The pairs ((x) r , E) constitute our 
estimate of the effective potential and are plotted in Fig.Q] 
for C = (full lines) and C = 0.01 (dashed lines). The 
agreement between the semi-classical estimate and numer- 
ical results is reasonable. 

The perturbation £ = 0.01 breaks the inversion sym- 
metry, therefore the right minima in g > g d regime in this 
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Figure 2 Average displacement and displacement fluctua- 
tions obtained by NRG compared to the semi-classical es- 
timate x m (dotted). 



case is lower in energy. In this regime, the molecule will 
predominantly reside near the right lead. Conversely, for 
g < <7d the potential is only slightly perturbed. 

3 NRG results Having established well the emergence 
of the double well effective potential we now check how it 
is reflected first in the static properties and then the dy- 
namical response of the system. To obtain these quanti- 
ties, we have performed the numerical simulations using 
the well known Wilson's numerical renormalization group 
1 26 27) (NRG) method. We restrict ourselves to the limit of 
zero temperature (T — > 0). The details of the calculations 
are given in Ref . (9) . 

3.1 Static quantities We begin by looking at the static 
quantities. The average displacement (x) for £ = van- 
ishes (as expected for an operator of odd parity under in- 
version in a state of well-defined parity). The fluctuations 
of displacement {Ax 2 ) 1 / 2 = ((x — (a;)) 2 ) 1 / 2 , shown in 
Fig.|2](full line) increase monotonically with g. The slope 
of (Ax 2 ) 1 / 2 is increased considerably at g ~ g d [or (~ 
g c ), see Ref. J9)], where the double well like effective po- 
tential is formed. This change of slope is driven by the in- 
creased hybridization in the odd-channel. 

For ( = 0.01 the absence of inversion symmetry is re- 
flected in the nonvanishing average displacement (dashed- 
dotted), which monotonically increases with increasing g. 
Therefore the fluctuations of displacement (dashed) in this 
case reach a maximum and then decrease with increasing 
9- 

For comparison, we plot also the position of the min- 
imum of the potential x m (dotted line) obtained from the 
semi-classical estimate Eq. (|6j, 



4^ 



(9) 



3.2 Phonon propagator Now we turn to the renor- 
malization of the phonon propagator by the electron-phonon 
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Figure 3 Displacement spectral functions for ( = 0. 
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Figure 4 The frequency of the soft mode peak as a function 
of g. The weight of the soft mode peak for £ = 0.01 and 
normalized to some arbitrary value (dotted). 



coupling. The dynamical information about oscillator is 
contained in the displacement Green's function. The dis- 
placement spectral function 



-Im <C x, x ^> u 



1 

7T 

i r°° 

-Im / (-i)([x(t),x(0)]}e iUJt dt (10) 
I" Jo 



is an odd function of u due to the hermiticity of x. Since 
A(ui) is in NRG evaluated for a finite system it consists 
of several rj-peaks of different weights. To obtain a smooth 
spectral function we have used the Gaussian broadening 
on the logarithmic scale [28], where the Dirac S function is 
broadened according to 



5{bJ - LO n ) 



1 



■ exp ■ 



log(w/w„) 



(11) 

and we used b = 0.3 in our calculations. 

In Fig. [3] we plot A(lu) for various g and £ = 0. The 
width of the high frequency peaks is overestimated due to 
the broadening procedure described above (for example, 
the width of the peak at co = Q for g = should vanish). 
We could use the Dyson equation [29 , 30 1 to obtain sharper 
peaks but we avoid this complication because on one hand 
there is no a priori guarantee that such a procedure gives 
more accurate results for large g and on the other hand in 
this work we are interested only in the position and not the 
width of the peaks. 

For intermediate g (starting at g ~ 0.5 for the param- 
eters used here) the vibrational mode begins to soften; the 
characteristic frequency of the oscillations is decreased. At 
still larger g > g^ two peaks emerge. The high frequency 
peak corresponds to the oscillations within each of the min- 
ima of the double-well potential, and the low-frequency 
peak (we denote its position by ujq) corresponds to the slow 
tunneling between the degenerate (or near-degenerate for 
( > 0) minima. 




Figure 5 Conductance, ( = 0.01. Shaded area (g > g c ) 
indicates the unphysical regime. 



The propagators for finite £ and £ = look alike, pro- 
vided g is small enough that ujq does not decrease below 
the frequency given by the energy difference oc £ between 
the minima of the two wells. 

For larger g the high frequency behaviour remains sim- 
ilar but as shown in Fig. fflwo which decreases for £ = 
exponentially (full line), for finite ( (dashed) saturates to 
the value oc ( of the energy difference between the min- 
ima. In this regime, the tunneling of the oscillator between 
the minima as characterized by the weight of the soft-mode 
peak in the phonon propagator (shown dotted) is suppressed. 

3.3 Conductance We show the conductance calcu- 
lated by NRG in Fig. [5] In the particle-hole symmetric 
point and at zero temperature the system is in the unitary 
limit, hence the conductance (in units of Go) is for g = 
near unity (reduced only due to small breaking of inver- 
sion symmetry for a value of order £ 2 ). For increasing g it 
decreases and at g ~ g^ if drops to zero. This occurs be- 
cause the molecular orbital is increasingly hybridized only 
to the right lead and the coupling to left lead V(l — gx) be- 
comes small. The point where the conductance is zero cor- 
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responds to the decoupling of the left lead, V (1 — gx) ~ 0, 
on the average. 

Surprisingly, for g larger still, the conductance increases 
again. This is due to the overlap with the negative ampli- 
tude V{ 1 — gx) < when states for which 1 ± gx < con- 
tribute considerably to the transport. As discussed in more 
detail in Ref. (9) this regime is an artifact of linearization 
and cannot be observed in the break-junction experiments. 
Here we remark only, that in the large-g regime the spin 
is screened by the odd channel. At the delimiting value 
of g (which we denote by g c ) the spin is simultaneously 
screened by the both channels and the resulting state is 
characterized by the non-Fermi liquid two-channel Kondo 
fixed point in the NRG flow. 

4 Conclusions In this report we considered a metal- 
lic break-junction bridged by a molecule. We were inter- 
ested in the electron transport through the break-junction, 
which is influenced by the modulation of tunneling be- 
tween the molecular orbital and the leads due to the os- 
cillations of the molecule with respect to the leads. 

We have shown that due to the electron-phonon cou- 
pling the harmonic potential confining the molecule to the 
center-of-inversion evolves to the double well effective po- 
tential. The change in the form of confining potential is re- 
flected in dynamical properties - e.g. a low-frequency peak 
emerges in the phonon propagator - as well as in the static 
properties - e.g. the fluctuations of displacement are in- 
creased. 

The emergence of the double well effective potential 
makes the system strongly susceptible to perturbations break- 
ing the inversion symmetry. When the frequency of the 
soft-mode decreases below the energy scale of such a per- 
turbation, the molecule is attracted to one of the leads. As 
a consequence, the conductance is suppressed. 

There are many points which we have not discussed in 
detail, including the breakdown of the linearization and the 
physics of the two-channel Kondo fixed point. The more 
comprehensive analysis of this interesting system will ap- 
pear elsewhere (9). 
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